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The time evolution of a set of 22Mq unstable charged stars that collapse is computed integrating 
the Einstein-Maxwell equations. The model simulate the collapse of an spherical star that had 
exhausted its nuclear fuel and have or acquires a net electric charge in its core while collapsing. 
When the charge to mass ratio is Q/VCM > 1 the star do not collapse and spreads. On the other 
hand, it is observed a different physical behavior with a charge to mass ratio 1 > Q/y/CM > 0.1. 
In this case, the collapsing matter forms a bubble enclosing a lower density core. We discuss an 
immediate astrophysical consequence of these results that is a more efficient neutrino trapping 
during the stellar collapse and an alternative mechanism for powerful supernova explosions. The 
outer space-time of the star is the Reissner-Nordstrom solution that match smoothly with our 
interior numerical solution, thus the collapsing models forms Reissner-Nordstrom black holes. 
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I. INTRODUCTION 


In this work we intend to study the effects of an 
electric field on the collapse of a massive star. We per¬ 
form direct relativistic simulations assuming spherical 
symmetry and integrating the Einstein-Maxwell equa¬ 
tions. We studied the collapse of several stars with 
different values of the total electric charge and de¬ 
termined the limits at which the collapse is avoided 
by the effect of the repulsive electric field. The elec¬ 
tric charge is carried by the particles that compose 
the collapsing fluid. The electromagnetic field and 
the internal energy of the gas contributes to the total 
mass-energy of the star, so it is not clear whether it is 
possible to overcharge a collapsing configuration. The 
analysis and conclusions drawn for collapsing charged 
shells cannot be directly extrapolated for this case, 
and a complete relativistic calculation performed in a 
self-consistent way is needed to know the outcome of 
the charged collapse. 

In particular, we found the formation of a bubble 
that was not predicted before. Although its theoreti¬ 
cal explanation is quite natural, the formation of the 
bubble depends on the initial conditions and its evo¬ 
lution is far from obvious due to the non-linearities of 
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the Einstein-Maxwell equations. Numerical solutions 
for the problem of charged stellar collapse were found 
only recently [l^, [d^ . 

It is commonly assumed that the stars are nearly 
neutrally charged due to a mechanism of selective ac¬ 
cretion of charges from the surrounding interstellar 
medium. Therefore, if a star is initially positively 
(negatively) charged, accretion of negative (positive) 
charges from the surrounding gas will tend to neu¬ 
tralize the total net charge. The same reasoning is 
applied to black holes [i^ . We observe that the effect 
of selective accretion were never fully calculated to our 
knowledge. On the other hand, an star can be electri¬ 
cally neutral while possessing a huge internal electric 
field [H]. 


The internal electric field of an astronomical ob- 
iect can be very large in some special scenarios [d^, 
1^. d^. In particular, the charge separation 

inside a spherical compact star, in hydrostatic equi¬ 
librium, can be very large when one of the plasma 
components is a degenerate gas while the other is a 
Maxwell-Boltzmann gas, i.e., like the gas of degener¬ 
ated electrons and the gas of nuclei in a white dwarf 


The calculations presented in this paper corre¬ 
sponds to the case in which the star posses a net elec¬ 
tric charge. However the results can be taken as an 
approximation to the more complicated problem of an 
star with total zero charge and with a non zero inter¬ 
nal electric field. We leave this point aside and for 
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clarity we will concentrate here on the dynamics of an 
star with a total positive charge. 

We ask whether an electrically charged star can 
collapse to form a charged black hole. Would the 
Coulomb repulsion avoid the collapse of the star ?, 
and moreover, are there physical differences with the 
physics of uncharged collapse ?. Moreover, from pure 
analytic analysis is not clear if it is possible to form 
an overcharged black hole (with Q > y/GM). The 
dynamics of a collapsing star could be quite different 
from the collapse of a charged shell onto an already 
formed Reissner-Nordstrom black hole, mainly due to 
backreaction effects. In the present paper, we will try 
to give the answer to these questions. 

Extremely charged black holes [Q/^/GM = 1) rep¬ 
resent an extreme limit in the context of the cosmic 
censorship hypothesis, since bodies with charge equal 
or higher than extremal are undressed by event hori¬ 
zons and constitutes naked singularities [1^ , [i^ . On 
the other hand, the formation of Reissner-Nordstrom 
(RN) black holes, is in connection with the third law 
of black hole thermodynamics i, da, di , because 
an extreme RN black hole has zero temperature. So, 
there is an interest on the formation of charged black 
holes from a pure theoretical point of view [49|| . 

In this paper we are not concerned with the mecha¬ 
nism that produce charge or an internal electric field. 
Neither magnetic fields nor rotation are taken into ac¬ 
count, although they could increase the magnitude of 
the electric fields. 

Without loss of generality, the density of charge was 
chosen proportional to the rest mass density. The 
interior solutions found can be matched smoothly with 
the Reissner-Nordstrom exact solution for the vacuum 
space-time dl- Therefore, in the cases in which the 
star collapses the result is the formation of a Reissner- 
Nordstrom black hole. 

The paper is organized as follows. In Section II, 
we discuss the problem of charged collapse analyzing 
the order of magnitude of the physical quantities. In 
the Section III are described the general relativistic 
equations that govern the dynamics of the stellar core 
collapse. In Section IV, are described the numerical 
techniques used to integrate the general relativistic 
equations and some caveats on its applications. Sec¬ 
tion V brings a discussion and interpretation of the 
numerical results. In the Subsection A, we show the 
calculation of the maximum mass formula for Newto¬ 
nian charged stars, and in Subsection B we calculate 
the optical depth for the neutrinos emitted during the 
collapse. We also discuss the implications for core 
collapse supernova. We end in Section VI presenting 
some final remarks. 


II. SOME ORDER OF MAGNITUDE 
ESTIMATES 

In this section we show some order of magnitude 
estimates, in order to put in a clear perspective the 
problem of a self-gravitating charged fluid sphere. The 
calculations in the present section are valid only in the 
non-relativistic regime. 

The formation of a charged star, is possible when 
the gravitational attraction overwhelms the electro¬ 
static repulsion on each single particle of the gas that 
is collapsing, i.e.. 


^grav ^ Select i 

or equivalently 

qQs >mMs , 

where Qs and Mg are the charge and mass of the star 
and q and m are the charge and mass of the particles. 
In this section we use a Newtonian approach, the full 
relativistic case could be different. 

Assume a star that has a total charge to mass ratio 
Qs/VGM s- According to Eq. ([2]) any particle with a 
charge to mass ratio q/VGrri can be added to the star 
if 

- (vfc) ’ 

where the equality gives the maximal charge to mass 
ratio 


( 1 ) 

( 2 ) 


Oi 



( 4 ) 


Eor example, a proton has a charge to mass ratio 
qp/'/Gnip ~ 10^®. Using the Eq. ([3]) the proton can 
be added to the charged star if 

<10"^®. (5) 


VGMs 


Thus, a star with higher charge to mass ratio can 
not be assembled from protons alone. On the other 
hand, if the infalling particles are dust particles with 
a larger charge to mass ratio, the maximum limit for 
the charge to mass ratio of the charged star can be 
much higher. 

An example of this is a self-gravitating ball of Eel 
nucleus, of mass Mhaii and charge Qbaii- As the charge 
to mass ratio for Eel is q/^/Gm ^ 10^®, the Eq. ([3]) 
for Fe gives 


0.2 


Qball ^ 10 “^® 
VG Mball 


( 6 ) 
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and we see that in this case a 2 ^ ai (greater than for 
protons alone, see Eq. ([5]) above). 

For “charged dust” formed by particles with 
ql^/Gm ^ 1, a self-gravitating sphere with mass 
Mdust and charge Qdust, can be in hydrostatic equi¬ 
librium if: 


Q!3 


Qdust 


VG Mdust 


1 . 


(7) 


U\- 

M- 


The mechanism produce a charge imbalance of 


~ 10 — charges per baryon. (10) 


This is several orders of magnitude higher than the 
needs for producing an extremal star. Therefore, 
charged old stars must be considered and studied as 
a theoretical possibility. 


So : as ^ 02 3> ai. Let us estimate the amount 
of charge needed to have an extremely charged 10 Mq 
star. In this case Qstar/VG Mstar = 1, so 


III. EQUATIONS IN CO-MOVING 
COORDINATES 


Qstar = V6.67 10-8 X 10 X 1.9 10^3 

= 4.9 X lO^o statCoulomb. (8) 

Therefore, it is needed to have an excess of 


We want to solve the Einstein-Maxwell equations 
(see Ref. [l3|l: 
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( 11 ) 


^ Qstar! % = 10"^° charges, 

in the star, where = 4.803 x 10“^° statCoulomb is 
the charge of one proton. There is roughly ~ 10^® 
baryons in a IOMq star. So, there must be one 
charged particle on ~ 10^® neutral ones in order to 
have a maximally charged compact object. This is 
compatible with the limits given above. Therefore, in 
the extremal case, there are a small charge imbalance 
equal to: 


dQ ~ 10 ^® charges per baryon, (9) 

in the star. A tiny amount of charge from a micro- 
scopic point of view, but with a huge total sum (see 

Ref. dl). 

It is possible to make the objection that the nu¬ 
cleus would disintegrate, or suffer nuclear hssion, be¬ 
fore assembling a charged star. However the energy 
locally available (in the center of mass of the particles) 
from the collapse is not enough to unbind a nucleus 
(bounded by the strong force). This was calculated 
for a compact charged neutron star. An exception is 
a nearly extremal compact object, in this case is ener¬ 
getically favorable for the nucleus to disintegrate [l^ • 
The conclusion is that we can approach to the forma¬ 
tion of a “nearly extreme” object, although it is not 
possible to reach the extremal value, the particles will 
fission first (see 0 for details). 

Studying theories with extra dimensions, Mosquera- 
Cuesta and co-workers found that particles lo¬ 
cated in the brane can leak out to the bulk space. It 
results that electrons can leak out more easily than 
baryons, producing a charge asymmetry that can be 
very large in very old stellar systems. This is suitable 
for type II supernova progenitors and neutron stars 


where is the Ricci tensor, is the metric ten¬ 
sor, R is the scalar curvature, and is the energy 
momentum-tensor, composed by two parts: 

T^^'' = T^'^ + ( 12 ) 

the perfect fluid energy-momentum tensor Tp'^ and 
the Maxwell energy-momentum tensor . The elec¬ 
tromagnetic energy-momentum tensor is given by: 

^ , (13) 


here is the Maxwell electromagnetic tensor, which 
can be written in terms of a potential A^: 

= A^.p - A^.p . (14) 


In the last equation and in the rest of the paper, we 
use semi-colon to denote covariant derivative. The 
perfect fluid energy-momentum tensor is: 


= {P + 8 c^) uV -Fg^^', 


(15) 


where 8 (? is the density of mass-energy and P is the 
pressure, given in [dyn/cm^], and vF is the 4-velocity 
of the observers co-moving with the fluid [1^, (il ]. 
0 . The electromagnetic tensor must satisfy the 
Maxwell equations [sj, [1^, 0: 


F^!^ = . (16) 

From Eq. d and using the Bianchi identities, we 
have: 


7^[a/3;7]=0. (17) 

We will use the common split of the energy: 

8 = p(\F (.!<?), (18) 
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where p is the rest mass density, and p e is the internal 
energy density [I^, [s^. 

The fluid must satisfy the energy-momentum con¬ 
servation equation: 


= 0 , (19) 

and the baryon conservation equation: 

{pu'')-,y = Q. ( 20 ) 

We will write the Einstein-Maxwell equations in co¬ 
ordinates and gauge choice appropriated for its nu¬ 
meric implementation [l3| . Considering a spheri¬ 
cal star, the line-element in co-moving coordinates is 
given by 

ds^ = a{t, — b{t, pL^dp'^ — R{t, p)‘^dfl‘^. (21) 

The 4-velocity of an observer co-moving with the fluid 
is u'' = [a“^, 0 , 0 , 0 ], and satisfies u'^Ui, = c^, where 
c is the speed of light. The coordinate p was gauged 
to be the rest mass enclosed by co-moving observers 
standing on spherical layers of the star. Each layer has 
its own constant value of p. In particular, observers 
at the surface of the star will be designated with a 
coordinate p = ps- 
We define [fiojl 


u = R^t/a . 


So, the equation of motion, obtained from the 
Einstein-Maxwell equations [l^ is 


u,t = -a 


9 r 

47ri?^ — 
w 



QQ,A 

47ri?4 J 


Gm 




47rG„ 

C2 c2 i?3 


, ( 22 ) 


where ri; = l-|-e/c^-|- Pj p<? is the relativistic specific 
enthalpy, p is the rest mass density, e is the inter¬ 
nal energy per unit mass, P is the pressure, G is the 
gravitational constant, Q{p) is the total charge inte¬ 
grated from the origin of coordinates up to the spher¬ 
ical layer with coordinate p, and m = m(t, p) is the 
total mass-energy defined below. We use the notation 
f,x = df/dx, with X = {p,t}, throughout the paper. 

We note that the Eq. resembles a Newtonian 
equation of motion plus relativistic corrections. In the 
Eq. (1^^ the last two terms are pure relativistic, and 
the full equation is equivalent to its Newtonian coun¬ 
terpart when c ^ oo [5l| . The factor T is a general¬ 
ization of the special-relativistic 7 factor for general 
relativity and is given by 



2mG 

Rc^ 


GQ^ 


(23) 


This factor also verifies the equality T = R fj,/b (see 
Refs, dl, [13, i). 

1 /2 

The equation for the metric component a = g^j. jc 
is obtained from the energy-momentum conservation 
equation = 0 , and is given by 0 


(qw)^^ _ 1 , p . Q Q,n 

aw wc^ \PJ fj. 47ri?^p 


(24) 


The other independent component Tq.^, = 0 gives 


„ = (25, 

which is identical to the non-relativistic adiabatic en¬ 
ergy conservation equation, and constitutes the first 
law of thermodynamics. We observe that the Eq. 
do not contain any electromagnetic term. This is due 
to the symmetry of the problem, although an electro¬ 
magnetic term arise in the total energy (see the Eq. 
Ell below). 

From the conservation of the number of baryons, we 

1/2 

obtain the equation for the metric component b = 

0 , 


b{t, p) 


1 

At: R? p 


(26) 


The Lagrangian coordinate p can be chosen to be 


nR' 

P = At: / pR^dR/T, (27) 

Jo 

this is the total rest mass enclosed by a sphere of cir¬ 
cumference 2tt R'. Thus, the collapsing ball is divided 
in layers of constant rest mass p and each co-moving 
observer is at rest in each of this layers [s^ - 

The equations for the total mass and the mass con¬ 
servation are, respectively d 3 : 


m{t, p) = At: 




RJ R^fj, dp + 


(28) 


1 


c 


2 



QQ,^ 

R 


dp; 


{pR^),t/pR^ = . (29) 

The charge 4-current is the product of a scalar 
electric charge density pch times = pchu"■ 

It can be shown that the charge and the rest mass 
are conserved in shells co-moving with the fluid [ 3 , 
0 . So, the charge increment between layers can be 
written as being proportional to dp, 

dQ = 4 TT Pch R^ dR/T , 


c 4 R 2 ■ 


(30) 
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where Pch oc p. 

We observe that Eq. dSOD is a non-linear equation, 
since E depends on the integral value of the charge 
Q{p) (see Eq. [531), so set of equations above must 
be solved self-consistently. 

The equations (I22l)-(l30l) above are solved numeri¬ 
cally with the following boundary conditions: 


The collapse05v2 uses a leapfrog method plus a 
predictor corrector step, and iterates with a Crank- 
Nicholson algorithm. The method is second order ac¬ 
curate in time and space. We use a numerical viscosity 
to resolve the strong shock waves formed [^ . 

A. Initial conditions and numerical caveats 


P = 0, 

at ^ = ps, 

vt, 

(31a) 

0=1, 

a.t p = ps, 

Vt, 

(31b) 

It = 0, 

at t = to) 

V/r, 

(31c) 

i? = 0, 

at p = 0, 

Vt, 

(31d) 


where ps is the mass coordinate at the surface of the 
star. The boundary condition (I31ap can be derived 
from the matching condition between the interior and 
the exterior solution [I^. Equation (|31bl) express our 
coordinate freedom to choose the time synchronized 
with a co-moving observer moving with the surface of 
the star. Eq. (l31cD means that the ball is initially 
at rest (at the initial time to)) and Eq. (I31dl) is the 
impenetrability condition at the origin of coordinates 
(this condition can be stated equivalently as m = 0 at 
/r = 0). In the next section we will present the initial 
conditions. 

In this work we choose 

Pch = constant x p , (32) 

for simplicity and without loss of generality. 


IV. NUMERICAL TECHNIQUES AND CODE 
SETUP 


The equations given in the section above are written 
in a form closely paralleling the equations of May and 
White [l^, [3, for the non-charged case [l^, p5j |. 
0 - 

We built the numerical code collapse05v2 to inte¬ 
grate the Eqs. (1551) - (1301) . 

The Eq. (1531) is re-written as 


aw 


aoWo exp 



QdQ \ 
pj 




(33) 


This equation is integrated numerically from the sur¬ 
face of the star toward the center, using the boundary 
condition agWs = 1, where and Wg are given at the 
surface of the star (with coordinate /i^), note that we 
already choose Us = I (see Eq. I31bl) . 


It is assumed a charge density proportional to the 
rest mass density through the star: pchlVCp ^ 
constant. From Eqs. dSZl) and dsop we see that 
Q{p)/VCp = constant, V (t, p). 

We assumed a polytropic equation of state, 
P = kp'^ , where 7 = 5/3 is the adiabatic coefficient 
[ 5 ^ . The star has initially a uniform mass density, 
p = 10^ g cm“^ and a uniform distribution of electric 
charge density, in all the models studied. The initial 
specific internal energy e is varied in the different sim¬ 
ulations. The initial setup imitates a massive star of 
M = 2 I.O 35 M 0 that had exhausted its nuclear fuel. 
Its fate depends on the internal energy and charge of 
the initial configuration. This is explored with our 
simulations. The initial radius of the star is Rq = 10® 
cm, in all the simulations. 

Also we consider that the electric field is always 
much lower than the Schwinger field limit for pair cre¬ 
ation (~ 10^® Volt/cm, in vacuo and in flat space-time 
[^1. For example, in the case of extremal charge 
Q/y/GM = 1, the total charge will be Q ~ 10^^ stat- 
Coulomb. In this case, the maximum electric field 
strength is ~ 10^^ stat-Volt/cm or ~ 3x 10^® Volt/cm. 
However, during the collapse the field can strengthen 
and take higher values, so quantum effects would be 
taken into account in a more realistic approach. 

We are not considering magnetic fields in the 
present simulations [s^. 

The simulations converge to the same solution when 
the resolution or the viscosity parameter are changed. 
The parameters of the simulation were chosen for ac¬ 
curacy and efficiency when running on a single pro¬ 
cessor. We performed the simulations using an Intel 
Pentium IV 2.6 Ghz processor (32 bits) and compiled 
the code with Lahey/Fujitsu v6.I [2l| for 32 bits; also 
we used an AMD Athlon EX 2.6 Ghz processor (64 
bits) and compiled with Absoft [I| for 64 bits, run¬ 
ning under Linux operating system. The code also 
runs under the Windows operating system and com¬ 
pile with the Fortran 90 Power Station. We didn’t find 
essential differences between the different runs with 
the exception of the speed of the simulation, which is 
faster with the AMD processor. In addition we set: a) 
the number of integration points: 1000; b) viscosity 
parameter: 3.3; c) initial time-step: 10“® sec. We also 
performed tests simulations with 300 and 1500 points 
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to be sure of the physical results, and tested the code 
varying the viscosity over a wide range of parameters. 
In the present simulations we tried to minimized the 
effects of the viscosity while keeping the ability of the 
code to resolve shock fronts. 

The percentage of numerical error in the Hamilto¬ 
nian constraint, expressed by the conservation of the 
total mass-energy is given by: 

Eh = 100 X |to — m'\lm, 

where m' = is the exact value of the total mass 

energy (constant in time) and m = m(t, /Xg) is the nu¬ 
merical value. The momentum constraint is “build 
in” the code and exactly conserved. In the present 
simulations: Eh < 0.1% when the AH forms, using 
1000 points, and it is Eh < 0.009 % with 1500 points, 
with no essential qualitative differences between the 
solutions. This values of Eh are for the worst case of 
a large binding energy and with very strong shocks, 
so it can be considered as the superior bound of the 
error. 

After crossing the AH the error grows to: ~ 0.4% 
(using 1000 points), and the simulation is stopped. 
Thus in every simulation we set a maximum error tol¬ 
erance of Eh ^ 0.4 %. We checked that using a higher 
number of points the errors are reduced. 

We emphasize that the numerical errors (mainly 
Eh, not the roundoff errors) and the intrinsic diffi¬ 
culty on the simulation is the formation of the black 
hole region, because near it several numerical opera¬ 
tors blows up. 


B. The apparent horizon 

For each layer of matter with Lagrangian coordinate 
/i it is possible to define an internal and an external 
horizon, 

»'±(t,/i) = ± 

(34) 

which can be interpreted as a generalization of the 
Schwarzschild radius for charged stars. During the 
simulation an apparent horizon forms when some 
layer of matter cross its external horizon, i.e. when: 

M- In addition to the forma¬ 
tion of the apparent horizon, in some cases a coor¬ 
dinate singularity develops, i.e.: 0 and 

- 1-00 (for some coordinate value /i' 0). 

Note that at the coordinate p, = 0 is located the 
physical timelike singularity of the vacuum Reissner- 
Nordstrom space-time. The simulation must be ended 
before any observer g can meet its internal horizon 


r_(/r). Hence, we can not say if the star collapses di¬ 
rectly to a singularity or if it pass through a wormhole 
to another asymptotically flat universe. In this paper, 
we leave this question open. 

The formation of an apparent horizon (AH) is a 
sufficient condition (although not necessary) for the 
formation of an event horizon and a black hole region 
[T^. Moreover, the interior solution matches with the 
exterior Reissner-Nordstrdm spacetime H- In this 
sense, we can say that the outcome of the simulations 
in which an AH forms is the formation of a “Reissner- 
Nordstrdm black hole”. Of course, if not all the matter 
collapse to the AH the result at the moment the sim¬ 
ulation is stopped will be an electric black hole plus 
a surrounding gas with a vacuum Reissner-Nordstrdm 
exterior spacetime. 

We observe that there must exist other formulations 
of the problem in which it is possible to get closer to 
the singularity while avoiding the AH. We will explore 
this subject in future works. 


V. DISCUSSIONS OF THE RESULTS 

We performed a set of simulations on the collapse 
of a 21 Mq star varying the initial conditions. We 
found that for certain cases a shell-like structure of 
higher mass density and charge is formed in the most 
external part of the star. This shell enlarge with time 
until it reaches the center of the star. 

The Table I summarize the results of the simula¬ 
tions performed. In each model studied are different 
amounts of total electric charge; initial internal energy 
or binding energy. The fate of each model depends on 
this conditions. 

The formation of the shell is sensible to the internal 
energy, or the binding energy, of the star (see Table 
I). When stronger is the bound of the star higher is 
the wall of the shell. 

In the Fig. o, we show the profiles of the mass 
density versus the Lagrangian mass g for the model 3 
of the Table I, in this case the shell forms only mildly. 
However, there is a strong shock wave formed. The 
Fig. © shows the velocity profiles, on which it is 
seen a strong shock propagating outwards. We ob¬ 
serve that although the shock is propagating outwards 
the sign of the fluid velocity is negative everywhere in¬ 
dicating that the star is always collapsing. There is 
an exception in which the shock velocity acquires a 
positive sign, although it lasts a very brief period of 
time. This is shown on the inset of the Fig. m-- the 
shell is reaching the coordinate origin, and an instant 
later it impinges the center. After that a shock wave 
is formed and start to propagate outwards. During a 
brief lapse of time the velocity is positive, although 
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the shock is not strong enough to stop the total col¬ 
lapse. The Figure ([2]) contains the snapshots of the 
metric component it can be seen that concor¬ 

dant with the shell impinging the center, gu acquires 
a value greater than one, indicating a blueshift with 
respect to the infalling matter. The inset of the Fig. 
([ 21 ) shows the details of the jumps produced on the 
value of g^l at the shell and at the shock positions, 
as well as the mentioned blueshift. 

Figure (|3]) contains the snapshots of the mass den¬ 
sity for the model 11 of the Table I. In this case there is 
an strong shell formed. On this figure we can see that 
the shell propagates inwards, to lower values of the La- 
grangian coordinate /r, until it impinges the center. At 
this moment, and place, it is formed an strong shock 
that starts its outward propagation. The inset box 
of the Figure (|4]) is a zoom of the dotted box, show¬ 
ing the details of the collapsing shell and the shock 
formed. The Fig. shows the verification that the 
collapse produce a black hole. In spherical symmetry 
the apparent horizon is formed simply when an ob¬ 
server cross his/her own Schwarzschild radius. The 
inset box of the Fig. ® shows the radius profile at 
the moment that the curve touch the Schwarzschild 
radius profile at a point. That point indicates an ob¬ 
server crossing the Schwarzchild radius and is the locus 
of the apparent horizon. After then the apparent hori¬ 
zon will evolve approaching the event horizon at in¬ 
finite coordinate time or after complete collapse (not 
calculated). 

Figure ®, shows the velocity profiles for model ff, 
there is an strong-shell and strong-shock in this case. 
However, the velocity has a negative sign and the mat¬ 
ter collapse directly to a black hole. Figure © shows 
the speed profile at the moment of time in which the 
shell impinges the center, and the formation and prop¬ 
agation of the shock wave. The inset box shows a 
close up of the details. There is a short lapse of time 
in which the shock has positive velocity. The Fig. ([S]) 
shows the metric coefficient g^l going to zero. In this 
case gtt never becomes greater than one, and the in¬ 
ner matter is always redshifted respect to observers 
standing at the surface or outside the star. 

Figure ([9]) shows the density snapshots for the col¬ 
lapse of an uncharged star (model 13). It is possible to 
compare this simulation with the former simulations, 
we see that for the charged collapse the density profiles 
are flatter all the way through the black hole forma¬ 
tion. Figure 113 shows the velocity profiles (compare 
it with the Figs. [2] and [6]). The Figure (fl^ shows 
the collapse of gl/'^ (compare it with the Figs. [2] and 
0). Figure m shows the profiles of the metric func¬ 
tion R, the bold line is the Schwarzschild radius at 
the moment the apparent horizon forms. The inset 


box shows the detail of the radius profile touching the 
Schwarzschild profile at one point, signaling the for¬ 
mation of an apparent horizon at the point indicated 
with the vertical arrow, at the time indicated on the 
legend box. The simulation is continued beyond that 
event, as seen in this figure. 

It is possible to compare the profile of the mass- 
energy function for the charged collapse with very 
strong shocks (model 3) Fig. (fT21) . with the case of 
uncharged collapse (model 13) Fig. iMl) . The inset 
box of the Fig. pA]) shows a zoom of the right end 
of the curves which indicates the level of energy con¬ 
servation in the simulation. All curves must reach the 
same point for perfect energy conservation. At the 
moment of the apparent horizon formation the mass- 
energy is conserved with a precision of 99.98 %, which 
is an excellent value taking into account that this is 
the worst case. After the apparent horizon forma¬ 
tion, the good conserving property is lost as in any 
relativistic code, and the simulation is ended with a 
0.48% error. We emphasize that this is the physical 
quantity we use to check the convergence properties of 
our code, since the momentum constraint is conserved 
exactly by the algorithm. As was said before the sim¬ 
ulations were performed using 1000 points, and the 
conservation can be better using a larger number of 
points. The Fig. Cl represents the energy conserva¬ 
tion for the case of uncharged collapse. The precision 
is roughly the same as in the former case, although 
using only 100 points. We checked that using 1000 
points the errors are lowered by a factor of 10 in this 
case, i.e.: 99.998% precision at the apparent horizon 
formation. The reason for better conservation of the 
energy constraint is that there are not very strong 
shocks like in the former simulations. 

The explanation of the shell formation can be grasp 
from the equation of motion (Eq. [22]). This effect 
can also happen in Newtonian physics, but its evo¬ 
lution in the strong field regime is highly non-linear 
and far from obvious. It seems difficult to make a 
complete analytic description of its evolution. How¬ 
ever, starting with a ball of constant rest mass den¬ 
sity and charge, it is easy to see that a shell must be 
formed. The relativistic term GQ^jc^ R^, in the Eq. 
(l22ll . produce a repulsion of the matter on the initial 
homogeneous sphere, which in turn produce a charge 
gradient QQ^fj^jAirR*. 

The term QQ AtkR^ in the equation of motion can 
be important in supporting the weight of the star (see 
Eq. 1221). Near the surface of the star, the pressure 
gradient term cancels the charge gradient term 
at a certain point we call /rg. Since at 
this point the two terms cancels out, all the matter 
outside the layer of coordinate /ig is almost in free 
fall, producing an accumulation of matter and giving 


Density plot 
model 3,Q/G’> = 0.484 


rise to the shell, while at coordinates with fj. < the 
matter is supported by a positive net gradient. We 
took care to check that this is not a numerical artifact, 
changing the boundary conditions at the surface, the 
amount of charge, and the initial conditions. 

Hence, it results that the point fig is roughly 
the locus in quo the term P^(/ro) and the term 
QQ,ij./‘iT^R‘^{n-o) sum zero, then the gas is compressed 
and a shell of matter forms. 

The non-linearity in the strong held regime comes 
from the fact that T is a function of the charge and the 
mass and multiplies the two gradient terms, enhancing 
the effect (see Eq. [22l) . 

From all the numerical experiments we performed, 
we observed that the shell formation effect arise more 
clearly when the initial density prohle is hat. For the 
case of charged neutron star collapse the ehect is neg¬ 
ligible [III- 

The point /xq shift to lower values with time, and 
this brings the shell to enlarge until /xq = 0, when 
it reaches the origin and rebounds forming the shock 
wave mentioned above. It is observed that the shock 
do not form in the uncharged case, using the same set 
of initial conditions. The density contrast between the 
shell and its interior is higher for greater total charge 
(or equivalently higher internal electric held and zero 
total charge). This is because the Coulomb repulsion 
inside the shell will be higher, and this gives a higher 
compression of the matter at the shell. The Coulomb 
repulsion term together with the pure relativistic term 
aGQ^/c^R^ in the Eq. ((2^ . are responsible for the 
positive mass density gradient near the coordinate ori¬ 
gin. The charged collapse is slightly delayed respect 
to the uncharged one, i.e.:~ 0.75 sec for model 1 re¬ 
spect to the uncharged case that lasts ^ 0.66 sec. At 
the end of the simulation, we observe that the density 
prohle looks globally hatter and the maximum density 
value is lower than in the uncharged case. 

The case 14 in Table I is special because it repre¬ 
sents stars with a charge greater or equal than the 
extremal value. We verihed that in this cases the star 
expands “forever”, i.e.: we follow its expansion un¬ 
til the density gives a machine underhow. In any of 
these cases we expect a re-collapse of the matter: the 
binding energy per nucleon tends to zero when the 
charge tends to the mass (see M), and is negative 
for a charge greater than mass. 


A. The maximum mass limit 

For charged stars there must exist a maximum mass 
limit like for uncharged compact stars. For white 
dwarfs the limit is known as Chandrasekhar mass 
limit, while for neutron stars is the Oppenheimer- 



FIG. 1: Snapshots of the mass density versus the mass 
coordinate /x, for the model 3 of the Table I, the shell forms 
mildly in this case, although the shock is very strong. 

Speed plot 



FIG. 2: Snapshots of the speed versus the Lagrangian 
mass, the inset shows the details of the shell impinging 
on the center and the shock starting its propagation. The 
speed is positive during a short lapse of time. 


Volkoff limit, etc. In this section, we calculate the 
mass limit for Newtonian charged compact stars. We 
assume, for simplicity, an equation of state dominated 
by electrons [^ . 

P/p = YekT/ruB + . 


(35) 
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Metric plot 
model 3, Q / = 0.484 



^ [solar masses] 


Density plot 
model 11, Q/G''" [1 = 0.290 



FIG. 3: Snapshots of the coefficient it is distinguished 

a sharp discontinuity at the shock and shell positions. The 

1 /2 

inset shows a zoom of the dotted box, the coefficient 
is greater than one at the instant of time the shock revert 
the sign of the velocity, indicating a blueshift respect to 
the infalling matter. The shock and the shell are indicated 
in this figure. 


FIG. 4: Snapshots of the mass density versus g, for the 
model 11 of the Table I. In this case there is an impor¬ 
tant imploding shell. The dotted box is zoomed on the 
inset, which shows the shell impinging at the center and 
the shock forming and propagating outwards. 


For a spherical configuration, the Newtonian hydro¬ 
static equilibrium gives: 

Pc/Pc = GM/R - QV(47ri?V) ■ (36) 


Simplifying the two equations, for T = 0 and 7 = 4/3 
we get (GM^ — Q‘^)/M'^^^ = ■ From this 

equations, we obtain the limiting mass Mmax- As¬ 
suming the charge is a fraction a of the mass, i.e.: 
Q = ay/GM, 




K 


G(l- 


a2) 


3/2 




5.83 


1 ? 


(l-a^) 




© ! 
(37) 


which reduces to the known Chandrasekhar mass for¬ 
mula when a = 0 (no charge), modulo some geometric 
factor depending on the density profile. 

For the extreme case a = 1, the configuration dis¬ 
perse to infinity. We checked that this also happens 
in the relativistic case (see model 14 of Table I). 


B. Implications for core collapse supernovae 

The quantity of energy deposition behind the 
formed shock wave is critical to produce a successful 
supernova explosion, like in the proposed mechanisms 
for core collapse supernovae explosions 0. El- 


In the cases we studied the shock wave is not fast 
enough to reach the surface of the star before the ap¬ 
parent horizon forms ahead of it, and the shock ends 
trapped by the black hole region. However, when in¬ 
cluding neutrino transport in our simulations the shell 
formation could have a dramatic effect. 

This provide a new mechanism for a successful core 
collapse supernova explosion. Although, the calcula¬ 
tions performed in the present work are far from a 
complete supernova calculation, it is simple to show 
that in the charged collapse there are a more efficient 
energy deposition to revive the supernova shock wave. 

As the density profile is flatter than in the un¬ 
charged case, the neutrino-sphere and gain radius for 
neutrino deposition must be different. In order to 
show this, we consider the change in the radius Ri,, of 
the neutrino sphere, which is defined by t{Ri,) = 2/3, 
where the optical depth of the uncharged matter for 
the neutrinos is 0 

t{R) = 1.5 X 10“ {N'^/QA) el i?-^. (38) 

According to the numerical results presented, we make 
the approximation p ~ constant^ for the charged case, 
in order to obtain an analytic expression for the opti¬ 
cal depth for neutrinos. In this case, the optical depth 
of the charged matter is 

TchiR) = 10-® (7VV6A) el pi2 (Rs - R). (39) 

It can be seen that Tch ^ t. Assuming typical val¬ 
ues 0: Cjy = 10 Mev and p ^ 2 x 10^°gcm“^, in the 




























































10 


Radius plot 
model 11.Q/G''"-^l = 0.29 



p, [solar masses] 


FIG. 5: The radius (LogioRadius) versus the Lagrangian 
mass, the bold line shows the Schwarzschild radius at 
the instant of the apparent horizon formation, the inset 
box shows the detail of an observer crossing his/her own 
Schwarzschild radius (an arrow is pointing the observer) 
and signaling the formation of an apparent horizon. It can 
be appreciated that the simulation continues far beyond 
the formation of the apparent horizon. 


uncharged case the position of the neutrino-sphere is 
at: Rn ^ 110 km, and in the charged case it 
is: Ri, ^ 845 km. This implies a more efficient neu¬ 
trino trapping during the charged collapse, and hence 
a more efficient energy deposition. 


VI. FINAL REMARKS 

We have performed relativistic numeric simulations 
for the collapse of spherical symmetric stars with a 
polytropic equation of state and possessing a uniform 
distribution of electric charge. We have not studied 
in this paper how the matter acquires a high internal 
electric field before or during the collapse process. 

In the present model we studied the essential fea¬ 
tures of the collapse of stars with a total electric 
charge. The simulations are an approximation to the 
more realistic problem of the temporal evolution of 
an star with an strong internal electric field and total 
charge zero. This model is under consideration by the 
authors. 

In the cases were the stars formed an apparent hori¬ 
zon it is unavoidable the formation of a Reissner- 
Nordstrom black hole. This is guaranteed by the sin¬ 
gularity theorems and by the Birkhoff theorem. 


Speed plot 

model 11,Q/G‘'V = 0.29 



FIG. 6: Snapshots of the speed versus p for model 11, it 
can be observed a jump in the velocity profile at the shell 
position dislocating through the center of the star. See 
Fig. 0. 


Speed plot: zoom mode 
model 1I.Q/G‘'V = 0.29 



FIG. 7: Detail of the speed’s profile showing the moment 
of time the shell impinges at the center of the star, later 
the shock is formed and starts its outward propagation, 
the inset shows a zoom of the dotted box. See Fig. ®. 


For low values of the charge to mass ratio 
q/Vgm < 1, no difference with the collapse of 
an uncharged star was found. The value of the to¬ 
tal charge that prevent the collapse of the star with 
any initial condition is given by a charge to mass 
ratio Q/VGM > 1. That means that stars with 
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TABLE I: 

Simulations of the collapse of a massive star varying the charge to 

mass ratio and the binding energy. 

Simulation 

charge to mass ratio 

Q/a/Gui 

[ergs/c^] 

Binding energy*" 

collapse ? 

shell ? 

1 

0.484 

9.61 X 10“'' 

-2.994 X lO""* 

yes 

yes 

2 

0.484 

1.92 X 10“® 

-2.992 X 10"^ 

yes 

no 

3 

0.484 

2.88 X 10"® 

-2.990 X 10~2 

yes 

no 

4 

0.484 

1.44 X 10“® 

-b2.498 X 10“^ 

no 

no 

5 

0.484 

1.35 X 10“® 

-1.764 X 10“® 

no 

no 

6 

0.145 

1.15 X 10"® 

-3.830 X 10“^ 

yes 

no 

7 

0.145 

9.61 X 10"* 

-3.833 X 10"^ 

yes 

no 

8 

0.145 

9.61 X 10“® 

-3.833 X 10“^ 

yes 

yes 

9 

0.290 

9.61 X 10"'^ 

-3.583 X 10"^ 

yes 

no 

10 

0.290 

9.61 X 10“® 

-3.584 X 10“^ 

yes 

yes 

11 

0.290 

9.61 X 10~“ 

-3.585 X 10"^ 

yes 

yes 

12 

0.850 

9.61 X 10“’^ 

-1.073 X 10“^ 

yes 

yes 

13 

0 

9.61 X 10"'^ 

-3.740 X 10"^ 

yes 

no 

14 

> 1.0 



no 

no 


^The energy density is initially distributed uniformly on the 
star in all the simulations. The equation of state is P = k 
with 7 = 5/3. The mass of the star in every simulation is 
M = 21.035 M©. 

^The binding energy of the star: (m — p.) in units of M© 
ergs. 


Metric plot 
model 11,Q/G‘'^H = 0.29 



(1 [solar masses] 


FIG. 8: Temporal snapshots for versus /r for model 
11 and its late collapse to zero. We must observe that the 
apparent horizon forms before gtt 0, which indicates 
the convergence towards the event horizon. 


Density plot 

q/g''V=o 



|t [solar masses] 


FIG. 9: Snapshots of the mass density versus /i for model 
13, an uncharged star. 


Q/\/GM > 1 spreads and do not collapse to form 
black holes nor stars in hydrostatic equilibrium. It 
is observed a dramatically different physical behavior 
when Q/^/GM >0.1. In this case, the collapsing mat¬ 
ter forms a shell-like structure, or bubble, surrounding 
an interior region of lower density and charge. The ef¬ 


fect is due to the competence between the Coulomb 
electrostatic repulsion and the attraction of gravity. 
This effect must occur in non-relativistic physics as 
well. The relativistic case is more interesting due to 
the non-linearity of the Einstein equations. The effect 
is more important when lower is the internal energy of 
the star (see Table I). In some of the experiments we 
observe a blueshift produced at the bouncing shock 









12 


Speed plot 

q/g''V=o 



FIG. 10: Snapshots of the speed profile, where u has 
the meaning of the rate of radius change per unit proper 
(comovel) time change. 


Metric plot 
Q/g'''-]i = 0 



FIG. 11: Temporal snapshots for versus for model 
13. 


wave (see Fig. H, because the strong shock wave ac¬ 
quires a positive velocity over a small lapse of time. 

For all the cases studied, the density prohle is glob¬ 
ally flatter than in the uncharged collapse. 

The optical depth for neutrinos is much higher in 
the case of the charged collapse, hence, the neutrino 
trapping must be more efflcient. In conclusion, if the 
internal electric field increase during the stellar col¬ 
lapse the formation of the shell must be taken seri- 


miQpr in GnnnrnmT'Ci cimniaCirinG 


Railitis plot 
Q/G‘”^l = 0 



g [solar masses] 


FIG. 12: The radius versus g: the bold line shows the 
Schwarzschild radius at the instant of time the apparent 
horizon is formed, the inset shows the detail of an observer, 
pointed with an arrow, crossing his/her own Schwarzschild 
radius, indicating the formation of an apparent horizon at 
this coordinate and at the time indicated. 


Newtonian charged stars, which clearly precludes the 
formation of naked singularities. The mass limit can 
not be naively extrapolated for the relativistic case, 
because charge and pressure regeneration effects can 
change the maximum charge to mass ratio. However, 
in the present work, we checked that it is not possible 
to form an star from matter with a charge to mass 
ratio greater or equal than one in agreement with our 
mass formula and with the cosmic censorship conjec¬ 
ture. 

A more complete simulation including neutrino 
transport, and other quantum effects, is being con¬ 
sidered by the authors. 


Acknowledgments 

The authors acknowledge the Brazilian agency 
FAPESP. This work was benefited from useful dis¬ 
cussions with Dr. Samuel R. Oliveira and Dr. Ed¬ 
uardo Gueron. We acknowledge Gian Machado de 
Castro for a critical reading of the manuscript. The 
simulations were performed at the Parallel Computer 
Lab., Department of Applied Mathematics, UNI¬ 
CAMP (Campinas State University). 































13 


Total mass-energy v.? Lagrangian mass 
model 3, Q / G''^ n = 0.484 



FIG. 13: Snapshots for the mass-energy versus /r for 
model 3. The numerical error in the simulation of the 
models can be measured by the failure of the curves to 
converge at the same point, this is indicated on the zoom 
boxes. After the apparent horizon formation the error is 
0.02 % for model 3 and model 13. The same error is found 
for model 13 using a lower number of radial points (see 
text for details). Compare with Fig. (I14II . 
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There are in the literature numerical and analyti¬ 
cal solutions of charged spherical stars in hydrostatic 
equilibrium (See for example Refs. [13 , [U, @l)- 
There are also numerical research on the collapse of 
charged scalar fields, with particular emphasis on the 
mass inflation effect near the Cauchy horizon [3^ (and 
references therein), and global anal ysis on the mass in¬ 
flation effect on the stellar collapse [3^ . The physics of 
charged thin shells collapse had been analyzed in 
and . The equations for stellar charged collapse had 
been obtained in the references (ah [2^, and specially 
for their numeric integration in |14|. 

The total charge of a star can be zero while its internal 
electric field is huge. An example of this are charge 
separation effects in compact stars, see Refs. [2^, [2^; 
charge separation during accretion of matter onto a 
star or onto a black hole, see Ref. [3^ and ■ electric 
field generation due to magneto-rotational effects [U]; 
electric field in pulsars [2i or not charge effects due 
to physics in extra-dimensions [T7| . 

Electromagnetic fields plays an important role in sev¬ 
eral models for the inner engine of gamma ray bursts. 
During the accretion of plasma onto a black hole the 
electromagnetic field can annihilate above a certain 
limit due to quantum effects, and produce a “dyado- 
sphere” of electron-positron pairs that in turn anni¬ 
hilates to produce a gamma ray burst [^ . 

This is the u^ component of the co-moving observer’s 
4-velocity in Schwarzschild coordinates, see Ref. M 
for details. 

In the hydrostatic limit we can take u^t = 0 in the 
Eq. (1221) to obtain the Tolman-Oppenheimer-Volkov 
equation for charged stars 0. 

The definition of the Lagrangian mass coordinate is 
independent of the presence of electromagnetic fields 

0- 

For 7 = constant shocks see [^, pag 128. 

A warning must be given: a magnetic field is induced 
in the collapsing charged sphere if perturbed slightly, 
and the spherical symmetry of the problem is broken. 
However, the strength of the induced magnetic field 
can be neglected in a first approximation to the full 
problem, as long as \v'^ jf? ^ 1, with v the non-radial 
component of the 3-velocity of the fluid. In the present 
case, V = 0. 

See Ref. (38|| for the uncharged case. 
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[56] From the simulations, we obtain the radius Ra = 954 


km, at which p = 2 x 10^° gem for the particular 
case Q/\/GM = 0.48. 



